The Symmetrical Cascode OTA
Design
Christian Enz
Swiss Federal Institute of Technology (EPFL), Lausanne, Switzerland

Schematic of the symmetrical OTA.
This notebook presents the design of the symmetrical cascode OTA which schematic is shown in the above figure. Note that all nMOS transistors have an odd number and all the pMOS transistors an even number which explains the numbering process. We will design the circuit with M1a-M1b in a separate well.
All the parameters required for the design are given below. They correspond to a generic $0.18\,\mu m$ process.
The OTA specifications are included in the file specs.py. For this design, they are focused on the gain-bandwidth product and dc gain without any constraints on the flicker noise or corner frequency.
We need to start with choosing the curent gain $A_3$. We have seen that having a current gain $A_3 > 1$ improves the transconductance but this is at the cost of a higher thermal noise excess factor. We therefore will set $A_1=A_2=A_3=1$
Since $A_1 = A_2 = A_3 = 1$, transistors M2a, M2b, M2c and M2d are all identical and so are M3a and M3b. We therefore have the following variables:
In this example there is no specification on the thermal noise. Therefore the transconductance $G_{m1}$ is set by the gain-bandwidth product according to \begin{equation*} GBW = \frac{G_{m1}}{2 \pi C_{out}}, \end{equation*} where $G_{m1}$ is the gate transconductance of M1a and M1b and $C_{out}$ the total output capacitance \begin{equation*} C_{out} = C_o + C_L \end{equation*} with $C_o$ the parasitic capacitance at the output node and $C_L$ the load capacitance.
In order to minimize the input-referred noise and offset, the input differential pair should be biased in weak inversion. The transconductance $G_{m1}$ in weak inversion is then given by \begin{equation*}\ G_{m1} = \frac{I_b}{n U_T}. \end{equation*} The bias current $I_b$ is the current flowing in each transistor M1a and M1b when the input differential voltage is zero. The bias current provided by M5b is therefore $2 I_b$. The bias current must satisfy the following inequality: \begin{equation*} I_b \geq 2\pi n_{0n} U_T C_{out} GBW_{min}. \end{equation*} Neglecting the output parasitic capacitance $C_o$ for now, the minimum bias current to achieve the target GBW is given by
The corresponding slew-rate is given by
which we will consider as sufficient. If this would not be the case we need to increase the current or move to an adaptive OTA.
To have some margin on the $GBW$ to account for the parasitic capacitance at the output $C_o$ that adds to $C_L$ but which we cannot estimate at this point, we set $I_b$ to
and the inversion coefficient $IC$ to
The transconductance can be calculated from the $G_m/I_D$ function as
which leads to the following $GBW$
which is slightly higher than the target specification offering some margin.
Knowing the drain current $I_{D1}$ and the inversion coefficient, we can calculate the $W/L$ aspect ratio for M1a-M1b:
Before finalizing the sizing of the differential pair, we first will size the current mirrors.
The drain voltage of M2a and M2b should be set as low as possible for a given maximum common mode input voltage still keeping M1a and M1b in saturation. For a maximum input common-mode voltage $V_{icmax}$ given by
The minimum source-to-gate voltage of M2a-M2b $V_{SG2}$ for this maximum input common mode voltage $V_{icmax}$ is given by \begin{equation*} V_{SG2} = V_{DD} - V_{icmax} + V_{GS1} - V_{DSsat1}. \end{equation*} The gate-to-source voltage $V_{GS1}$ and the saturation voltage $V_{DSsat1}$ of M1a-M1b are given by
We set $V_{SG2}$ to
This corresponds to an inversion coefficient given by
The saturation voltage of M2a-M2d is then given by
The specific current transconductance and W/L are then given by
We will now have to make sure that the non-dominant pole $f_{p2}$ and $f_{p4}$ at nodes 2 and 4 are sufficiently higher than the GBW to insure the desired phase margin. The non-dominant pole is given \begin{equation*} \omega_{p2} = \frac{G_{m2}}{C_2}, \end{equation*} where $C_2$ is given by \begin{equation*} C_2 = 2(C_{GS2} + C_{GB2}) \end{equation*} Assuming M2 is in strong inversion and saturation, we have \begin{equation*} C_{GS2} \cong \frac{2}{3}\,W_2\,L_2\,C_{ox} + W_2 \cdot C_{GSOp} \end{equation*} and \begin{equation*} C_{GB2} \cong \frac{n_{0p}-1}{3\,n_{0p}}\,W_2\,L_2\,C_{ox} + W_2 \cdot C_{GBOp}. \end{equation*}
The capacitance at node 2 scales with $W_2$ and $L_2$ according to \begin{equation*} C_2 = W_2\,L_2 \cdot C_{WL} + W_2 \cdot C_W, \end{equation*} with \begin{align*} C_{WL} &= 2\,C_{ox} \left(\frac{2}{3} +\frac{n_{0p}-1}{3\,n_{0p}}\right),\\ C_W &= C_{GSop} + G_{GBop}. \end{align*}
Since the $W/L$ has already been set by the transconductance and the current, we can derive $W_2$ and $L_2$ for achieving a given capacitance $C_2$ according to \begin{align*} W_2 &= \frac{-C_W \cdot W_2/L_2 + \sqrt{W_2/L_2} \cdot \sqrt{4\,C_2\,C_{WL}+C_W^2 \cdot W_2/L_2}}{2\,C_{WL}},\\ L_2 &= \frac{W_2}{W_2/L_2}. \end{align*}
Setting the non-dominant pole $f_{p2}$ to 10 times the $GBW$, we get
We choose
We can check
Since $A_1=A_3=1$, M2a, M2b, M2c and M2d are all identical and hence
To size M3a and M3b we need to
Let's start maximizing the output swing assuming that we bias the cascode transistor M7 in weak inversion for minimum saturation voltage and maximum $G_m/I_D$. We set the inversion coefficient of M7 to
The saturation voltage of M7 is then given by
Let's say that we want an output swing of 1V this means that we have 0.8V to share between the positive and negative saturation voltages which we decide to split equally. This means 0.4V for the saturation voltages of M7 and M3b which leaves a saturation voltage for M3b:
which gives an inversion coefficient for M3b:
which we round to
We now can check whether with the chosen inversion coefficient for M3b, M2c remains in saturation. The $V_{GS3}$ is given by
Which is OK.
Knowing the current and the $IC$ we can deduce the $W/L$ as
Similarly to what was done for transistor M2a-M2c, we need to make sure that the non-dominant pole $f_{p3}$ at node 3 remains much higher than the $GBW$. Setting $f_{p3}/GBW=10$ and proceeding in a similar way we get
We see that $W_3$ is smaller than the minimum width. If we don't want to increase $W_3\,L_3$ and hence $C_3$, we need to reduce the inversion coefficient $IC_3$. We can find the $IC$ such that $W_3=W_{min}$ for the given $f_{p3}$ in the following way
The inversion coefficient of M3a-M3b is therefore reduced from 28 to 16 in order to maintain a reasonable parasitic capacitance at node 3 to make sure $f_{p3}$ is higher than the $GBW$.
M3b is identical to M3a
The cascode transistors are sized according to the desired DC gain given by \begin{equation*} A_{dc} = \frac{A_3 \cdot G_{m1}}{G_o} \end{equation*} with \begin{equation*} G_o = \frac{G_{ds3b} \cdot G_{ds7}}{G_{ms7}} + \frac{G_{ds2d} \cdot G_{ds4}}{G_{ms4}} \end{equation*} The DC gain is given by
and the output conductance is then given by
We will split the output conductance equally between the NMOS and PMOS cascodes. We already have set the inversion coefficient of M7
Knowing the current we can deduce $I_{spec7}$ and $W_7/L_7$ as:
Knowing the current and the $IC_7$ we can get $G_{ms7}$ as
We already have chosen the length of M3a and M3b
The output conductance of M3b can then be estimated as
We can then deduce the output conductance of M7 as
from which we deduce the length of M7
which smaller than the minimum length. We choose the minimum length
The bias voltage $V_{b2}$ has to be set such that M3b remains in saturation. The minimum bias voltage $V_{b2}$ is given by
For convenience we choose $V_{b2}=V_{DD}/2$.
M2d has already been sized according to:
The saturation voltage of M2d is then given by
Similarly for M4 we set the same inversion coefficient as M7
from which we get its saturation voltage
The output swing is then limited to the positive side by
Knowing the current and the $IC$ we can deduce $I_{spec4}$ and $W_4/L_4$ as
The source transconductance of M4 is given by
The output conductance of M2d can be estimated as
We can then deduce the output conductance of M4 as
from which we can calculate the length of M4 as
$L_4$ is smaller than $L_{min}$. We choiose $L_4 = L_{min}$
From which we get the width of M4
We finally choose
The bias voltage $V_{b1}$ has to be set such that M8 remains in saturation. The maximum bias voltage $V_{b1}$ is given by
We can set it to $V_{DD}/2$
We can check the dc gain
which is within specs.
Finally we need to size M5a-M5b. The saturation voltage of M5b is set by the minimum input common mode voltage
from which we can deduce the maximum saturation voltage of M5b
The corresponding $IC$ is given by
We choose
From $I_b$ and $IC_5$ we get $I_{spec5}$ and $W_5/L_5$ as
We need to set $W_5$ to the minimum width
We can round $L_5$ to
Transistor M5a is identical to M5b
We haven't finalized the sizing of M1a-M1b. We already have the $W/L$ ratio but still need to set $W_1$ and $L_1$. We don't have any specifications on the noise and particularly the corner frequency, while the white noise is already set by the transconductance $G_{m1}$.
Since the current mirrors current gains are all unity $A_1=A_2=A_3=1$, the transistors M2a, M2b, M2c and M2d are all identical. In this case the contributions of the current mirrors to the input-referred flicker noise relative to the contribution of the differential pair is captured by the $\eta_{fl}$ parameter recalled below \begin{equation*} \eta_{fl} = 1+2\,\frac{\rho_p}{\rho_n} \cdot \left(\frac{G_{m2}}{G_{m1}}\right)^2\,\frac{W_1 L_1}{W_2 L_2} + \left(\frac{G_{m3}}{G_{m1}}\right)^2 \cdot \frac{W_1 L_1}{W_3 L_3}. \end{equation*} The ratio of the flicker noise parameter of pMOS transistors to that of nMOS transistors is given by
which is rather large. The square of the transconductance ratio $(G_{m1}/G_{m2})^2$ is given by
which should compensate for this large $\rho_p/\rho_n$ ratio. We can try to make the contribution of the input differential pair equal to the contribution of the pMOS current mirrors
We choose
The contributions to the input-referred flicker noise of the pMOS current mirrors and nMOS current mirror relative to the contribution of the differential pair are given by
The gamma noise factor and the input-referred white noise are given by
This sets the corner frequency to
The specifications are recalled below.
The bias information are summarized below.
The transistor sizes are summarized below.
We can plot the magnitude and phase of the open-loop gain.
We now can estimate the input-referred noise PSD.
We see that the OTA thermal noise excess factor is only slightly larger than that of the differential pair. This is due to the rather large $G_{m1}/G_{m2}$ and $G_{m1}/G_{m3}$ ratios.
We can now compute the input-referred flicker noise and the corner frequency.
We see that the current mirrors contribute about the same than the differential pair.
We can plot the input-referred noise
The variance of the input-referred offset for $A_1=A_2=A_3=1$ is given by \begin{equation*} \sigma_{V_{os}}^2 = \sigma_{V_{T1}}^2 \cdot (1 + \xi_{V_T}) + \left(\frac{I_b}{G_{m1}}\right)^2 \cdot \sigma_{\beta 1}^2 \cdot (1 + \xi_{\beta}) \end{equation*} where \begin{align*} \xi_{V_T} &= 2\left(\frac{G_{m2}}{A_1\,G_{m1}}\right)^2 \cdot \left(\frac{A_{VTp}}{A_{VTn}}\right)^2 \cdot \frac{W_1 L_1}{W_2 L_2} + \left(\frac{G_{m3}}{G_{m1}}\right)^2 \cdot \frac{W_1 L_1}{W_3 L_3},\\ \xi_{\beta} &= 2\left(\frac{A_{\beta p}}{A_{\beta n}}\right)^2 \cdot \frac{W_1 L_1}{W_2 L_2} + \frac{1}{A_1^2} \cdot \frac{W_1 L_1}{W_3 L_3}. \end{align*} with \begin{align*} \sigma_{V_{T1}}^2 &= \frac{A_{VTn}^2}{W_1 L_1},\\ \sigma_{V_{T2}}^2 &= \frac{A_{VTp}^2}{W_2 L_2},\\ \sigma_{V_{T3}}^2 &= \frac{A_{VTn}^2}{W_3 L_3}, \end{align*} and \begin{align*} \sigma_{\beta 1}^2 &= \frac{A_{\beta n}^2}{W_1 L_1},\\ \sigma_{\beta 2}^2 &= \frac{A_{\beta p}^2}{W_2 L_2},\\ \sigma_{\beta 3}^2 &= \frac{A_{\beta n}^2}{W_3 L_3}. \end{align*}
We see that the $\beta$-mismatch is negligible and that the input-referred offset voltage is dominated by the contribution of the $V_T$-mismatch from the differential pair.
The total current consumption, ignoring the current drawn by M5a, is $I_{tot} = (2+A_1+A_3) \cdot I_b$. Assuming the input differential pait M1a-M1b is biased in weak inversion, the bias current is directly related to the gain-bandwidth product $GBW$ according to \begin{equation*} I_b \cong n U_T \cdot \frac{C_L}{A_3} \cdot GBW. \end{equation*} The total current consumption can then be written as \begin{equation*} I_{tot} \cong \frac{2+A_1+A_3}{A_3} \cdot n U_T \cdot C_L \cdot GBW. \end{equation*} The bias current and the total current are proportional to the gain-bandwidth product $GBW$ and the load capacitance $C_L$. In the case $A_1=A_2=A_3=1$, $I_{tot} \cong 4 \cdot n U_T \cdot C_L \cdot GBW$.
This notebook presented the design of the symmetrical OTA based on the design equations derived in the Analysis Notebook. The design will be checked by simulation in the Verification Notebook.